Hall Effect and Resistivity in High-Tc Superconductors: 
The Conserving Approximation 
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The Hall coefficient, _Rh, of high-Tc cuprates in the normal state shows the striking non- Fermi 
liquid behavior: _Rh follows a Curie- Weiss type temperature dependence, and |J?h| 2> i/\ne\ sit low 
temperatures in the under-doped compounds. Moreover, _Rh is positive for hole-doped compounds 
and is negative for electron-doped ones, although each of them has a similar hole-like Fermi surface. 
In this paper, we give the explanation of this long-standing problem from the standpoint of the nearly 
antiferromagnetic (AF) Fermi liquid. We consider seriously the vertex corrections for the current 
which are indispensable to satisfy the conservation laws, which are violated within the conventional 
Boltzmann transport approximation. The obtained total current Jk takes an enhanced value and 
is no more perpendicular to the Fermi surface due to the strong AF fluctuations. By virtue of this 
mechanism, the anomalous behavior of Rh in high-Tc cuprates is neutrally explained. We find that 
both the temperature and the (electron, or hole) doping dependences of Ru in high-Tc cuprates are 
reproduced well by numerical calculations based on the fluctuation-exchange (FLEX) approximation, 
applied to the single-band Hubbard model. We also discuss the temperature dependence of Rn in 
other nearly AF metals, e.g., V2O3, k-BEDT-TTF organic superconductors, and heavy fermion 
systems close to the AF phase boundary. 

PACS number(s); 72.10.Bg, 74.72.-h, 74.25.Fy 
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I. INTRODUCTION 

In the normal state of high-Tc superconductors 
(HTSC's), various quantities deviate from the conven- 
tional Fermi liquid behaviors. [Q These non- Fermi liquid 
features have been studied intensively both theoretically 
and experimentally because they have close relation to 
the mechanism of superconductivity. For example, the 
electrical resistivity (p) and the longitudinal NMR relax- 
ation rate (l/Ti) in HTSC show universally the behav- 
iors p (X T, 1/Ti cx T" for a wide range of temperatures. 

These are quite different from the conventional Fermi 
liquid behaviors, p cx T^, 1/Ti cx T. 

In HTSC's, the Hall coefficient (i?H) also shows an in- 
teresting non-Fermi liquid behavior: It shows a drastic 
temperature dependence although the Fermi surface (FS) 
in HTSC is non-degenerate and its shape is simple. At 
high temperatures (~ lOOOK), Ru takes a nearly con- 
stant value, and its doping dependence is very small. Its 
value is close to the one estimated by the LDA band cal- 
culation, R!^^'^. [|j The doping dependence of i?H is also 
very small there. 

On the other hand, as the temperature decreases, Ru 
begins to show a Curie- Weiss type temperature depen- 
dence, and its maximum value is a few times larger than 
jjband ^]-^g optimum-doping. This enhancement of 
i?H further increases in the under-doped region. In the 
hole-doped compounds, e.g., YBa2Cu3 07_5 (YBCO) or 
La2-5Sr5Cu04 (LSCO), dRn/dT < is observed and Rh 
is positive for a wide range of temperatures. [^-|9j On the 
other hand, dRB_/dT > is realized in the electron-doped 
compounds, e.g., Nd2-5Ce5Cu04 (NCCO), so the sign of 



i?H changes to negative at a low temperature although 
its FS is hole-like. [|jl|Jll| Figure shows a summary of 
experimental results of LSCO and NCCO in the under- 
doped region, where an approximate electron- hole sym- 
metry is realized. In both compounds, |i?H| increases 
near the half-filling. 
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FIG. 1. Temperature dependence of Ru in LSCO 
(hole-doping) and NCCO (electron-doping) in the paramag- 
netic state. Note that l/\ne\ ~ 1.5 x lO^^cm^C in HTSC's. 
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The Hall effect is one of the unsolved problems in 
HTSC. Its unusual features mentioned above are sum- 
marized as follows: (i) The Curie- Weiss type behavior of 
i?H in a quite wide range of temperatures, (ii) The en- 
hancement of i?H iu the under-doped region, (iii) Rn < 
in the electron-doped compounds. 

Nowadays, various non- Fermi liquid phenomena of 
HTSC have been explained by using different types of 
spin-fluctuation theories, e.g., the SCR theory, fl^] the 
spin-fluctuation model, Q and the FLEX theoryTl p^pSl 
They can explain a reasonable Tc of the dx2-y2 supercon- 
ductivity. They can also explain the pseudo gap forma- 
tion in the density of states, the shadow band 
formation, and the collective modes emerging below 

So far, various attempts have been made on the Hall 
effect in HTSC. Some of them are in the frame- 

work of the spin fluctuation model, by using the Boltz- 
mann transport approximation (or one-loop approxima- 
tion), p^p^ However, the Boltzmann approximation 
can not explain the magnitude of |i?H|(3> 1/ne) in the 
under-doped region. Moreover, it predicts i?H > for 
both YBCO and NCCO because they have similar hole- 
like FS's. This result contradicts with the experiments, 
shown in Fig. |l|. Therefore, the behavior of i?H noted as 
(i)-(iii) above have casted some suspicion on the validity 
of the nearly AF Fermi liquid description for the HTSC. 

In this paper, we study the Hall effect of HTSC's based 
on a conserving approximation. |25y2(]| ] We use the ex- 
pressions for the conductivity and the Hall conductivity 
derived from the Kubo formula. |p7|-^0| Then, we study 
the vertex corrections for the current according to the 
conserving approximation, and find out that they show a 
critical behavior as a natural consequence of the strong 
backward scatterings by the AF fluctuations. By virtue 
of this fact, we succeed in explaining the overall features 
of i?H, noted as (i)-(iii) above, without assuming a non- 
Fermi liquid ground state. We also find that the conven- 
tional Boltzmann approximation, where the conservation 
laws are violated, cannot reproduce any of (i)-(iii). 

We analyze the extended Hubbard model as an effec- 
tive model for HTSC's. We use the fiuctuation exchange 
(FLEX) approximation to calculate the Green function 
and the self-energy. (l^Jl^ It is a kind of self-consistent 
perturbation theory with respect to U, and it has advan- 
tages for handling large spin fluctuations. 

Phenomenologically, the spin propagator in HTSC's is 
expressed for small q and w as follows: |13,0,E3[ 



temperature as usual. In general, the following relations 
are satisfled for T > T* experimentally: pW 



XQ 



i + e2(q-Q)2 + iwM 



(1) 



where Q is the antiferromagnetic (AF) wavevector, and ^ 
is the AF correlation length. Experimentally, follows a 
Curie- Weiss type temperature dependence at higher tem- 
peratures. It ceases to increase at Tc in the over-doped 
region, or at the characteristic temperature T*(> Tc) in 
the under-doped region. We call T* the pseudo spin-gap 



(2) 
(3) 



where Q, ao, ai and 012 are constants. |32| ] The coef- 
ficient ao increases rapidly in the under-doped region, 
and reaches O(IOO) at T* nearby the half-fllhng. 
(We put the unit-cell length a — 1.) And a2 decreases 
rather moderately in the under-doped region. The re- 
lation ujsf >, T (wsf <, T) is satisfied in the over-doped 
(under-doped) region. 

The typical spin-fiuctuation theories reproduce the ex- 
perimental relations (|) and (|) for T > T*. ^ More- 
over, the approximate relations p cx £,^T^, 1/Ti fxf^T 
are derived in the nearly AF Fermi liquid. ||l^ , |23| , p^^ 
So, the experimental non-Fermi liquid behaviors of p and 
1/Ti are naturally explained by the spin-fluctuation the- 
ories. One may also expect that the anomalous behavior 
of i?H in HTSC comes from the T-dependence of In 
this paper, we find that i?H oc is realized through the 
vertex corrections for the current. A similar study based 
on the phenomenological AF spin-fluctuation model is 
reported in another paper, 

The contents of this paper are as follows: In §11, we in- 
troduce the single-band Hubbard model with some sets of 
parameters corresponding to YBCO, LSCO and NCCO. 
In §111, we review the general formulation for a^^ and 
Uxy/H based on the Fermi liquid theory, and rewrite 
Uxyl H into a simpler form. In §IV, the vertex correc- 
tions to the current is studied by using the conserving 
approximation. We find that only the Maki-Thompson 
term is dominant. In §V, we solve the Bethe-Salpeter 
equation for the total current Jk analytically, and the re- 
lation i?H oc is derived. In §VI, numerical results for 
p and i?H obtained by the FLEX theory are presented, 
which are consistent with the experimental behaviors in 
HTSC's. Finally, in §VII, the Hall effect in heavy fermion 
systems are discussed briefly. 

The readers who are mainly interested in the numerical 
calculation of i?H can proceed to §VI.B for the flrst read- 
ing, where the sufficient set of equations for the numerical 
calculations for Uxx aud Uxy/H are explained shortly. 



II. MODEL HAMILTONIAN 

In this paper, we treat the following extended Hubbard 
model with ([/, fo, ^i, ^2): 

H = Y1 «kcLck<. + XI 4+qT4'-qiCk'iCkT, (4) 

k(T kk'q 

= 2tQ[cos{kx) + cos{ky)) + Ati cos{kx) cos{ky) 
+2t2{cos{2kx) + cos{2ky)) , (5) 

where cj^.^ is the creation operator of an electron with 
momentum k and spin cr, and U is the on-site Coulomb 
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repulsion. We represent the filling of the electrons by n, 
and n = 1 corresponds to the half-filling. 
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Taking the results by the LDA band calculation into 
account, |^5|-^ we choose the following set of parame- 
ters: (I) YBCO (hole-doping), NCCO (electron-doping): 
to = -1, ti = 1/6, t2 = -1/5. i|,|6|,|9| (II) LSCO (hole- 
doping): to = -1, h = 1/10, t2 = -1/10. Figure 
H shows the Fermi surfaces (FS's) for U = together with 
those for finite U calculated by the FLEX approximation 
at T = 0.02. In the case of (I), the spectrum at (7r,0) 
is below the chemical potential /i at least for n > 0.6, 
and the FS is hole-like everywhere. On the other hand, 
in the case of (II) the spectrum at (tt, 0) is above /i for 
n < 0.77 at f7 = 0, and for n < 0.85 at C/ = 6. These 
characters of the FS's coincide qualitatively with those 
by the LDA band calculations ]3^-|3^ or by the angle 
resolved photoemission (ARPES) experiments. [[40[-j42[ 
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FIG. 2. Fermi surface of (a)YBCO, (b)NCCO, and 
(c)LSCO. In (a) and (b), to = -1, ti = 1/6, and t2 = -1/5. 
In (c), to = -1, ti = 1/10, and t2 = -1/10. 

Here, we summarize the formalism of the FLEX the- 
ory which will be used throughout this paper. The Dyson 
equation is written as 



{Gk(e„)} ^ = ie„ -I- - - Ek(e„). 



(6) 



(-71,-71) 



(IC.-7I) 



The self-energy I]k(eTi) given by the FLEX approxima- 
tion is 

Sk(e„) = T^Gk-q(e„ - Wi) • ^q(^i)' (7) 

^q(^0 = Qxq(c^/) + ix'^iu^l) - X° (C^O) + U, (8) 

xt,iu^i)^xl-{l~Uxl{coOy\ (9) 

Xq(^/) = X°-{l + t/x°(c^i)}"\ (10) 

X°(c^/) - -T^ Gq+k(^/ + e„)Gk(e„), (11) 



where e„ = {2n + l)7rT and oji — 21ttT, respectively. The 
self-energy is shown by Fig. We solve the equations 
(p|)-(pl]) self-consistently, choosing fi so as to keep the fill- 
ing constant, n = TJ2u n Gk(en) ■ e""^"''', where 6 = 0+. 
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FIG. 3. The self-energy of the FLEX theory. The full line 
and the wavy line represent G{k — q) and V{q), respectively. 
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n=1.15 



(7C.-7I) 



In a Fermi liquid, the real-frequency Green function in 
the vicinity of a; '-^ and |k| ~ fcp is represented as 



Gk(t^) = Zk/ (cj + /i - ek - i^kAk), 



(12) 



where Zk is the renormalization factor given by Zk = 
(1 — 9ReEk(ci;)/9u;)~^, ek is the quasi-particle spec- 
trum given by the solution of {ReGk(w)}^^ ~ 0, and 
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Ak = -ImEk(w + iS) > 0. The density of states (DOS) 
is given by 



Pk(w) 



—lmGii{uj 

TT 



iS). 



(13) 



In the case of ZkAk <C T, Pki^) — Z]^^ ■ 6(0; + fi — ek). 

The FLEX approximation is suitable for the analy- 
sis of the nearly AF Fermi liquid. It has been applied 
to the square lattice Hubbard model by many authors. 
@,[6|-[l8|j2^g| Though it is an approximation, imagi- 
nary time Green function obtained by the FLEX agrees 
with the results by the QMC simulations very well for a 
moderate U. Recently, it has also been applied to the 
superconducting ladder compound, Sri4_a;Caa;Cu2404i, 
and the organic superconducting k-(BEDT-TTF) 
compounds. 

We comment on the anisotropy of Ak on the FS, 
which becomes larger as the AF fluctuations grows at 
low temperatures. Ak takes a large value around the 
crossing points with the magnetic Brillouin zone (MBZ)- 
boundary, which we call hot spots as often referred to 
in literatures. And Ak becomes small at the 

points where the distance from the MBZ-boundary is the 
largest, which are called cold spots, (see Fig. |[) These 
cold spots play the major role for p and i?H- We study 
this subject in §VI in detail. 
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FIG. 4. The hot spots and the cold spots in YBCO and 
NCCO, respectively. 

Finally, we discuss the validity and the limitation of the 
FLEX theory on HTSC's. In the FLEX approximation, 
eqs. (H) and ^ are satisfied well, and the coefficient 
in eq. (|^) increases rapidly as n approaches to the half- 
filling (71 = 1). Moreover, the relation oj^f > T) in 
the under-doped (over-doped) region is satisfied qualita- 
tively as shown in Table I, which is consistent with exper- 
iments. However, the FLEX approximation can not re- 
produce the experimentally observed pseudo spin-gap be- 
haviors for T* > T > Tc, where ^ ceases to increase and 
l/wsf begins to decrease as the temperature decreases. 
It would also be inapplicable near the Mott-insulating 
state, i.e., 0.9 >, 1.1. In this paper we perform nu- 



merical studies for n < 0.9 or n > 1.1, where the FLEX 
approximation gives reasonable results. 

Table I : The value of Wst for n — 0.90(under-doped), 
n = 0.85(nearly optimum), and n — 0.80(over-doped) 
obtained by the FLEX approximation at T = 0.02. 





n = 0.90 


n = 0.85 


n = 0.80 


YBCO (C/ = 8) 


0.018 


0.034 


0.046 


LSCO (C/ = 6) 


0.013 


0.019 


0.024 



III. FORMALISM OF CONDUCTIVITY IN THE 
FERMI LIQUID THEORY 

In this section, we review the transport theory. By the 
Kubo formula, the conductivity is given by 



<7^ll^ — ^ ^k^^k'i/ 

»1/T 



ImXkcr,k'cr'(^ + 



kcr.k'cr' (i^ri) - 



UJ = 



X 



(Tr {ck<TWck<T(T)4'<x'CkV'}), 



(14) 



(15) 



where Wk^(fc) — de^/dk^ and ti>„ — 2TiTn is the even 
Matsubara frequency, and e(> 0) is the absolute value 
of the charge of an electron. In the absence of the mag- 
netic field, the analytic continuation from -ftTko-^kv (iw„) 
to i^Tkcr.k'cr' (fjJ + 1(5) has been performed by Eliashberg. 
[p7| According to him, the conductivity carried by the 
quasiparticles are given by 



5/ 

de 



1 



(16) 



where /(e) = 1/(1 + e^^-^'>''^). In eq. ^ we have 
done the energy integration by assuming the relation 
ZkAk ^ T, which is not always satisfied in HTSC as 
shown in §VI, however. 

In (|l|), Uka; and Jka; are given by 



d 



;e;^-fReSk(cc> = 0)) 



E 

k' 



^rkk'(0,e)|Gk'(e)|'- Jk': 
47ri 



(17) 
(18) 



where 7kk' (e, e') is the irreducible four point vertex, 
which plays an important role to treat the umklapp pro- 
cesses of conduction electrons. The total current Jk is 
given by the solution of the Bethe-Salpeter (BS) equation 
(p^), which is shown by Fig. |[ We note that %.\i'{e,e') 

is represented as 7^2°'' (^j ^'s w = 0) in the Eliashberg's pa- 
per. 13 
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FIG. 5. The BS equation for Jk^. 

The Hall coefBcient i?H under a weak magnetic field 
along z-axis H is give by 



<JxylH 



^ xx^ yy 



(19) 



The analytic continuation for the normal Hall conduc- 
tivity axy due to the quasiparticle contribution in the 
presence of the magnetic field H has been performed by 
Kohno and Yamada, |^ or Fukuyama et ai, in the 
gauge invariant manner. According to them, in case there 
is the four-fold symmetry of the system, 



k 

dk,. 



dl 
de 



^M-(k)- 



(20) 



where exyz = —^yxz = 1, reflecting the Onsagaer's reci- 
procity theorem. Equation ( pO| ) means axy/H oc (Ak)~^, 
whereas Gxx oc (Ak)~^ by eq. ([l^). Thus, i?H =const. 
in the conventional Fermi liquid at low temperatures. 

The expression ( po| ) can be rewritten into a simpler 
form, assuming that there is only the symmetry with re- 
spect to the origin, (see eq.(3.21) of ref. p^.) 



de 



2k 



(Ak)2 

J-kx ■ {Sz X i/k) V • Jky - {x ^ y) 



As(k) = Axy{k) + Ayxik) 



= \Vk\ 



T ^ T 

ok\\ 



(21) 



where is the component of k along the vector e\\ (k) = 
[cz X ?7k)/|'Ck|, and is tangential to the FS at k because 
Wk is perpendicular to the FS. In eq. (|2|), 0,/(k) is the 
angle between Jk and the x-axis, except for an arbitrary 
constant. Contrary to A^y(k), j4s(k) introduced in eq. 
( pl| ) is a scalar variable, i.e., independent of the choice 
of coordinates. As a result, axy/H is also independent of 
the choice of coordinates, if only the reflection symmetry 
exists. This property of axy/H has been proved so far 
only by the Boltzmann transport theory. [[46|-[48| 

By using the relation de^^/dk^ — Zk^k/^, eq. ( pT] ) be- 
comes 



axy/ H 



dk\\ I Jk 



FS 



|2 ( dem 

\ dk\\ 



(Ak)^ 



(22) 



at sufficiently low temperatures. In this line-integration, 
k-point moves anti-clockwise along the FS. 

Finally, we discuss on the Boltzmann transport the- 
ory: The conductivity in the magnetic field H is given 
by cr^,. = eX;k(-^//t^ek)i'kM*i^('^k) for E \\ e^, where 
*i.(wk) = (1 - eA-^Wk X i?) • V) • (eAj^^Wki.) up to the 
first order of H within the relaxation time approxima- 
tion. i pTj As a result, the conductivity in this approxi- 
mation, (7°^, is given by ( [T6| ) by replacing Jk with Uk. In 
the same way, the Hall conductivity within the relaxation 
time approximation is given by 



aljH^-- i dfcyKkl 

4 JfS 



dg.(k) 

dk\\ 



1 



(Ak)2^ 



(23) 



where 0^(k) is the angle between Wk and the a;-axis. 
Thus, the sign of axy/H is determined by the sign of 
dO^ (k) / dk\\ , which is nothing but the curvature of the 
FS at k. 11^] 

In the later sections, we calculate Jk^ by solving the 
BS eq. (|l^). In the nearly AF Fermi liquid, we find that 
Jk is no more perpendicular to the FS, so dO j / dk\\ and 
d6^{k.) / dk\\ at the same k can be quite different, even in 
sign. This is the reason why the Boltzmann approxima- 
tion fails to reproduce the anomalous behavior of i?H in 
HTSC's. 



IV. VERTEX CORRECTIONS FROM Tk,k'(e,e') 

In this section, we study the vertex corrections for the 
current, which is essential for the transport phenomena. 
The self-energy in the FLEX theory, which is given by eq. 
(^ , is also obtained by the functional derivative of <&flex 
as Ek(e) = (5$FLEx/<5Gk(e), where $flex is given by the 
closed skeleton diagrams made of Gk(e) and C/, with a 
factor 1/n for [/"-diagrams. The existence of $flex, 
which is depicted in Fig. ^, means that the FLEX theory 
is classified as a conserving approximation whose frame- 
work was constructed by Baym and Kadanoff, j2^ and 
Baym. In the conserving approximation, the particle- 
hole transport function L is given as the solution of the 
BS equation, where the irreducible particle-hole vertex 
rkk'(e,e') = <5Sk(e)/i5Gk'(e') is used as the kernel. Then, 
the L obtained in this way satisfies various conservation 
laws automatically. This is the reason why we call it the 
conserving approximation. 
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FIG. 6. Each full line represents the dressed Green function 
G\i(uo) and each broken line represents U . 
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Significance of the conserving approximation for the 
correlation functions is well recognized in various situ- 
ations. Conductivity is one typical quantity. Within 
the conserving approximation, Yamada and Yosida show 
that the conductivity given by (^) diverges at finite 
temperatures in the absence of the umklapp processes, 
reflecting the momentum conservation law. Their work 
also shows that the vertex correction from 7kk' (e, e') in 
( p^ ) , which is neglected within the Boltzmann theory, is 
indispensable to treat the umklapp scattering processes 
of electrons self-consistently. 

The irreducible particle-hole vertex Fkk' within 
the FLEX theory is shown in Fig. |^. 



Fkk'(e„,e„';c^/)=F('^)+FW+F(^), 



(24) 



where e„ — {2n + 1)ttT is the odd Matsubara frequency, 
and uji = 21ttT is the external frequency. We put the 
external momentum zero. F^^^-F^'^) are given by 



xGk+q(e„ +t^/')G'k'+q(en' 
rkk'(*^n,en';W/) = -T^ Wq(t^/',CJ;' + UJl) 

xGk+q(e,i + t^i' + t^O^k'-qlen' ~ ^^v 



(25) 



(26) 



(27) 



where Vk(en) is given by (||), and we have introduced 
Wk(wi,a;;') as 



W,,{louuv) = {Uxliu^i) + 1) {UxlM + 1) 

+ \u^ {Ux'k{u^i)-l){Ux'kM-l) 
-U^. (28) 

These three irreducible vertices F^"\ F*-^' and F'-'^-' are 
sufficient for the conserving approximation. In litera- 
tures, (a) process is called Maki- Thompson (MT) term, 
and (b) and (c) are called Aslamazov-Larkin (AL) terms, 
respectively. 
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FIG. 7. The irreducible four-point vertexes, which are suf- 
ficient for the conserving approximation. In (a), the wavy 
line represents V{k — k'). In (b) and (c), the two wavy hnes 
represent W{q) and the g-summation should be taken. 



In order to solve the BS equation (|18|), we have 
to obtain the functional form of the irreducible vertex 
^k'(e, e') in (P^. For this purpose, we perform the an- 
alytic continuation of F*^''^ F^*") and F^'^^ with respect to 
e„ and e„' in Appendix A. Thus, we obtain the vertex 
corrections for the current by replacing 7kk'(0,e) in eq. 

(|l|) with r^^,''\Q,e) given by eqs. (0)-(^. 



At first, the contribution coming from T-^]J is given by 



xpk'(e') 



Ak'(e') 



(29) 



where we put e = 0, and we have used the relation 
|Gk(e)P = 7rpk(e)/Ak(e). This vertex corrections play 
an important role in the singular behavior of the Hall 
coefficient in HTSC's, which will be discussed in §V. 



Next, we consider the correction terms come from 77 



(6) 



and T^j^./. Approximately, they are given by 

= I E / (^th^ - th^) f d^w, 

1,' n 



kk' 



(CO) 



X th 



2T 



th^j Pk+q(0)pk'+q(0) 



Ak'(O)"'*'''' 



(30) 



where Wq{uj) is introduced by eq. (|A6|). Equation 
( ^ ) is derived by neglecting the w-dependences of pi^iuj) 
and Ak(ti-'), which will be allowed qualitatively because 
only the regions |a;|, |e| <, minjT, o^sf | are important in 
the w, e-integrations in eq. (|30|). The variable change 
k' — k' is performed for AJ^. Strictly speaking, how- 
ever, AJ^{uj) and AJ^{uj) are not equal. 

Now, we show that AJ^^^^ is negligible in the case of 
^ I and Q = (7r,7r). In this case, the leading contri- 
butions in the q-summation in eq. (poh come only from 
q ~ Q. In the k'-summation of eq. pQ), there is a can- 
cellation between the contributions from k' and — k' if 
we put q = Q, because k' -I- Q and k' — Q are the same 
in the momentum space. As a result, AJ^'^ is expected 
to be negligible compared with AJ^^. This statement 
becomes rigorous in the case of —^ oo. In Appendix B, 
we show this cancellation in the two AL terms explicitly 
by the numerical calculations. This is one of the main 
results of this paper. 



V. ANALYSIS FOR THE BETHE-SALPETER 
(BS) EQUATION FOR Jk^ 

The aim of this section is to give the qualitative un- 
derstanding of the mechanism for the temperature depen- 
dence of i?H in HTSC's. We try an analytical approach 
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to solve the BS equation ( |T^ ) for Jk^, by neglecting the 
AL-terms. To simplify the discussion, we assume that 
Q — (tt, tt) and the MBZ-boundary lies on (±7r, 0) and 
(0, ±7r), which is realized in YBCO experimentally. 

For a qualitative discussion, we use the phenomenolog- 
ical expression for x^^{lo), eq. (j^), and neglect other terms 
in the definition of V^(w) in eq. (||). We introduce the 
function H{x) = l/x - '2.i){x + 1) + 2V'(x + 1/2), where 
^(x) is the di-Gamma function. Then, the imaginary 
part of the self-energy is given by 



(-It.IC) 



Ak = V I — [cth— - th— 



Imyq(e + \b) ■ Pk-q(e) 
(31) 
(32) 



\2ttTJ ~ Wq(Wq + 7rT/2) 



where Wq = LOgi + Wsf^^(q — Q)^. [g3[ We see 



that 

-1 



if (a;q/27rT) takes a large value only when |Q — q| <, ^ 
and it is negligibly small elsewhere. The approximate 
form of eq. (^ij) with eq. (^) is also obtained by extrap- 
olating the results of both limits, cjgf ^ T and cogf ^ T, 
where the e-integration in eq (31) becomes easier. 



Next, we examine the vertex correction for the current, 
eq. ( p9|) . We stress that ImVk(a;) appearing in eqs. ( |29| ) 
and (|31|) are same, which is ensured in the conserving 
approximation. We can show that 



sr Tjf Pk-q(O) , 

— ;- > , XQt^sf • H — ^ • Jk 



q 



V27rT/ Aw-r 



(33) 



Comparing eq. ( ^3|) with eq. (|3T|), and noticing that 
H{u!q/2TTT) is negligibly small for |Q — q| >, 1/C, we get 

Z\Jk W ( Jq )|q_k'|<l/J 

« Jk'-(cos(0,7(q)-0,7(k')))|q_k'|<i/^, (34) 



where k, k' and q are on the FS. Here, we have intro- 
duced k' on the FS so that {k'^,k'y) — (ky, k^) ■ sign(— fc^- 
/cj,), as shown in Fig. ^. We see the relation k — k' w Q 
is satisfied in the momentum space. Here we assume 
IQ — (k— k')| <, ^"^ around the cold spot because it seems 
to be satisfied in the present numerical calculation by the 
FLEX theory. Thus, we obtain a simplified BS equation. 



-'^k = + ak • ^k' 



(35) 



where ak ~ (1 — c/^^) < 1 and c '-^ 0{1) is a constant, 
ak takes the maximum value around hot spots. 



(K,JI) 



It' 




l(k 

\ , MBZ 


, ' k-k' / 


k-k' \ 


k'^\ 




YBCO 



X 



(-It-It) 



(Il.-It) 



FIG. 8. The relation between k and k'. Both points locate 
on the FS. 



Then, eq. (35) is easily solved as 
1 



l-«2 



(wk + ak • ""k') 



(36) 



Equation \Q>m means that Jk is not parallel to Uk. For 



example, (i) at k^ = ky, Jk — Wk/(l + ak) 



^iJk is satis- 



fied, (ii) Near the MBZ-boimdary, Jk « (Cv2c)(t/k+Vk')) 
which is nearly parallel or perpendicular to Q. (iii) On 
the BZ-boundary, Jk || Vk and Jk ~ Vk is realized be- 
cause the contribution from k'-point cancels out with 
that from —k'-point approximately, due to the fact that 
|k — k'l = |k + k'l w IQI in the momentum space. 
Thus, we should put ak = in eq. ( |3^ ) on the BZ- 
boundary. These behaviors of Jk together with Vk are 
shown schematically in Fig. H, which is confirmed by the 
numerical calculation in §IV B. Physically, this peculiar 
behavior of Jk comes from the multiple backward scat- 
tering between k and k' caused by the AF fluctuations. 



(0,71) 




(0,0) (7t,0) 
FIG. 9. Schematic behaviors of Jk and vu,. Contrary to vu,, 
Jk is not perpendicular to the FS. For example, {d6j/dk\\) < 
on the XY-axis, and {d6j/dk\\) > on the BZ-boundary. 
(see eq. (^.) 
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Apparently, Jk shows the divergent behavior near the 
AF phase boundary where ak ~ 1 is reahzed. Now 
we stress the importance of the conservation approxima- 
tion to get the correct ak- For instance, we always get 
oo if we replace |Gk(e)|2 with |G^(e)|2 in eq. ph, 



ctk 



which leads to the divergence of Jk. 

First, we consider the conductivity a^x- According to 
( |l6| ) , (T^a; is given by the averaged value of Wk^ • J^x over 
the FS, that is. 



{\Vkx\'^ - akIfkxWkyI} 



(37) 



We see that eq. (|37|)sa |wkxP/(l + ak) is realized around 
the cold spots in YBCO and LSCO, where jwka;! ~ l^kyl is 
satisfied. This means that the conductivity axx is smaller 
than that given by the Boltzmann approximation, due to 
the vertex corrections for the current given by 7kk'(^j ^')- 
This is confirmed by the numerical calculations in §VI. 

Next, we discuss the axy By using eq. (|3^), dJ]^x/dkf^ 
is given by 



dk. 



' Jkx 



Vky 



2ak/3i 



l-al 



1 



f dvwx 
V dkfj 



T ttk 



J]i.x 

dk,. 



(38) 



where /3k. = da^/dk^. Hereafter, the =F in equations is 
equal to sign{—k,j;-ky). By using eq. (p3|), which gives 
the definition of ak, then /3k, ^ is given as 



/?k,p~-e'(Q-(k-k')), 



(39) 



<, 1/^ is satisfied. We see that /3k, ^ 
^ is closer to the MBZ-boundary 



when IQ - (k - k' 
is positive when fc + e • e 
than k is (e^ is a unit vector along ^-direction, and e > 
is a small constant.), and vice versa. 

As a result, ^s(k) introduced in eq. (Elh is given by 



^ AJk) = — - 



Pk 



I'k X -TT-Vk 

afc|| 



=F/3k,|| 



(40) 



where the momentum fc|| is tangential to the FS, and is 
along the anti-clockwise direction. Thus, a^y/H is en- 
hanced by the factor, 1/(1 — a^) cx or /?k,|| oc 
contrary to the case of a^x- The first term of ( ^ ) is 
proportional to the contribution given by the Boltzmann 
transport theory, whose sign is determined by the curva- 
ture of the FS. It takes a larger value inside of the MBZ. 
(see fig. H.) On the other hand, the second term of pO| ) 
is negative inside of the MBZ, and is positive outside of 
it. This term is dominant outside of the MBZ because 
l/9k,||| 3> 1 and \v^^ — |wkP is satisfied there. We 

notice that /3k, || = on the MBZ-boundary and on the 
XY-axis. 

The obtained results in HTSC's with the strong AF 
fluctuations are summarized as follows (see Fig. ||.) : 



• The portion of the FS inside of the MBZ gives 
rise a positive contribution to Rn. In other wards, 
^(k) cx {dej{k)/dk\i) < inside of the MBZ. 

• The outside part of the MBZ gives rise a pos- 
itive contribution to i?H. In other wards, 
{de,j{k)/dk\\) > there. 



In the above ^j(k) is introduced in eq. (^2|). This change 
of the sign of i?H never occurs within the Boltzmann ap- 
proximation because {ddy(k) / dk\\) < everywhere. 

The Hall coefficient will be determined by the region 
near the cold spots because of the factor (Ak)~^ in eq. 
(p^). As shown in Fig. ^, the cold spots locate inside 
(outside) of the MBZ in the case of YBCO (NCCO). 
Thus, we understand the reason why Rn > in the hole 
doped systems, and why the sign of i?H changes in the 
electron-doped systems. 

In conclusion, the axy/H is proportional to cx 1/T 
both in the hole-doped systems and in the electron-doped 
systems. We find that the T-dependence of axx and 
axy/H for a system with the strong AF fluctuations are 



C^xx ^°/Ak, 

axy/H ^ ±e/Al 

i?H - ±e - ±XQ, 



(41) 



where +(— ) is for hole-doping (electron-doping) case. 
The factor comes from the vertex corrections for the 
current introduced in this paper, which does not ap- 
pear within the Boltzmann approximation. We will con- 
firm this analysis by numerical calculations based on the 
FLEX theory in the next section. 



VI. NUMERICAL RESULTS 



A. One-Particle Properties and Magnetic Properties 
obtained by the FLEX approximation 

Here, we show electronic properties obtained by the 
FLEX approximation. In this section, we use U = 8|to| 
for YBCO in numerical calculations, considering that 
the band width W is 8|to|. On the other hand, we 
use U = 6|io| for LSCO, and U = 5.5\to\ for NCCO, 
to reduce the Stoner factor in the FLEX calculation, 
as — maxq{?7 • Xq(0)}. We have checked the numeri- 
cal results do not depend on U qualitatively. 

In this section, we put \to\ — 1. Then, T — 0.1 will 
corresponds to ~ 600K because |to| 0.5eV in the LDA 
band calculation. In the present calculation, 4096 k- 
point meshes and 256 Matsubara frequencies are used. 
By solving the linearized Eliashberg equations, we obtain 
Tc « 0.02 (- 120K) for YBCO and LSCO at n = 0.85, 
which is close to the T^'s reported by the previous stud- 
ies based on the FLEX approximation. ||l5|,|l8|,|olE|] Also, 
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Tc « 0.01 for NCCO at n = 1.15. We find that the sym- 
metry of the superconducting state is dj.2_y2-\ike in all 
cases. 

Figure |lO| shows the temperature dependence of the 
FS's for YBCO, NCCO and LSCO, respectively. They 
are determined by the relation + Rel]k(0) = H- At 
low temperatures, the curvature of the interacting FS 
becomes smaller, which is more prominent in YBCO 
and LSCO. The deformation of the FS originates from 
Rel]k(0), whose analytical expression at zero temper- 
ature is given by ref. [ p4[ within the spin-fluctuation 
model. According to it, the sign of ReEk(O) is equal to 
sign(^ — Ek-q) approximately, which moves the FS to- 
wards the MBZ-boundary. As the temperature increases, 
the FS becomes closer to the non-interacting one because 
^ decreases. This temperature dependence of the FS 
should make |i?H| smaller as the temperature decreases, 
because its curvature around the cold spots decreases 
then. We examine this effect numerically in the next 
subsection. 

Moreover, in the FLEX approximation, the flat-band 
structure (i.e., extended saddle point) is created around 
the van-Hove singularity points, (±7r, 0) and (0, ±7r), be- 
cause of the renormalization effect by 1/zk ^ 10. |^] 
This is also the origin of the sensitive temperature depen- 
dence of the FS in YBCO and LSCO shown in Fig. |l^. 
This flat-band structure is actually observed by ARPES 
experiments. |Q 
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FIG. 10. The temperature dependence of the interacting 
FS's for various compounds. It is small and negligible in the 
case of (7 = 0. 

Next, we consider the q-dependence of the static mag- 
netic susceptibility Xq(0), given by eq. (p^. Because 
it does not contain the vertex corrections required in 
the conservation approximation, it gives a slightly over- 
estimated value in the under-doped region. |1§|| In real- 
ity, the observed Xq(0) by the neutron diffraction exper- 
iments can not be expressed by the simple spin-fermion 
model in eq. p): For YBCO, Xq(0) shows a peak around 
q « (7r,7r). | pO[ On the other hand, it is incommen- 
surate for La2_5Sr5Cu04, and shows a peak around 
q « ((1 - <5)7r, tt), (tt, (1 - S)t:) for 0.18 >S> 0.05. H 

Figure ^ shows the calculated (g^;, (7j;)-dependence of 
Xq(0) for YBCO (n = 0.90, T = 0.02), NCCO {n = 1.20, 
T = 0.02) and LSCO {n = 0.85, T = 0.06), respec- 
tively. We see that Xq(0) is commensurate for YBCO 
and NCCO, which is also consistent with neutron diffrac- 
tion experiments. In the case of LSCO, Xq(0) shows 
an incommensurate structure at low temperatures. At 
n — 0.85, the peaks locate at q = (0.837r, tt), (tt, 0.837r) 
at T = 0.02, which is consistent with experiments, and 
it becomes commensurate for T > 0.08. In conclusion, 
main characters of Xq(0) for each compounds are repro- 
duced well by the FLEX calculation with appropriate set 
of parameters, (to t2,U). 




qx (7t,0) 
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Figure |l^ shows the temperature dependence of 
maxq{xq(0)} of YBCO, NCCO and LSCO for differ- 
ent fiUing numbers. These plots are nothing but the 
T-dependence of Various experimental works on 

HTSC's by neutron diffraction or by NMR confirm that 
follows the Curie- Weiss law qualitatively for T > T*, 
and its Curie constant increases as n 1. [ p]j3^ As 
shown in Fig. |l^, the FLEX approximation reproduce 
both the temperature and the doping dependence of 
in HTSC's for T < 0.1. However, the calculation of Xq(0) 
including the vertex corrections will be required for more 
detailed studies. 



FIG. 12. The temperature dependence of maxq{x(q, 0)} 
for various filling numbers. All of them follow the Curie- Weiss 
law, which is a universal feature of HTSC's. In YBCO and 
NCCO, Xq(0) takes the maximum value at q = Q. We note 
that xq(0) oc 

Finally, we discuss the k-dependence of Ak — 
— Iml]k(0 + iS) on the FS. Here, we define 

A(fc||)= J dk^A^-pkiO)/ y"dfc±Pk(0), (42) 

where fc|| and k± are the momentum parallel and per- 
pendicular to the FS, respectively. A(fc||) is an averaged 
value of Ak over the fc^-direction on the FS, which has a 
finite width at finite temperatures. Figure O shows the 
fc|| -dependence of A(fc||) over the 1/8-part of the FS, as 
shown in Fig. In each cases, the relation ZkAk ^ T 
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is realized around the cold spots because 1/zk ^ 10 is 
satisfied. 
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FIG. 13. The fcy-dependence of A(fc||) at various tempera- 
tures. The T-dependence of A(fc||) at the cold spot and the 
hot spot are quite different. 



BZ (t.It) 




FIG. 14. The path of fc|| in the case of YBCO. 

For YBCO and LSCO, A(A:||) takes maximum value 
not on the hot spot shown in Fig. but on the BZ- 
boundary, because of the influence of the van-Hove sin- 
gularity around (tt, 0). As a result, the spectral weight at 
the Fermi energy is reduced around (7r,0), which is con- 
sistent with ARPES experiments. And A(A:||) takes min- 
imum at the cold spot. On the other hand, for NCCO, 
the hot spot locates on the XY-axis, and the cold spot 
locates on the BZ-boundary. It is quite important that 
the position of the cold spots changes across n « 1 by 
using the FLEX approximation, which causes the change 
of sign of i?H as shown in §V. To verify this result ex- 
perimentally, the ARPES measurements for NCCO are 
desired. 

As shown in Fig. [l3[ A(A:||) at the cold spots decreases 
in proportion to T qualitatively in all cases. To under- 
stand its behavior, we calculate the following quantity: 

(Ak)Fs = ^AkPk(0) 

k 

XQUJsi ■ W aimx^(^ + i<5) 



A 2^ ' 



4 ^ UJq{ujq + TrT/2) 

« XqTC^ • f 1 - (1 + 7rT/2c^3f) 



-1/2 



^=0 

(43) 



We notice that the q-dependence of ^ImXq(w + iS)\i^^Q 
in the FLEX approximation is moderate. Here we as- 
sume IQ — (k — k')| < even if k and k' = ±{ky, k^) 
are on the cold spots. Then, (Ak)FS is expected to reflect 
Acoid because pk(0) takes smaller values around the hot 
spots. As a result. 



Acold oc 



oc T 
cx T 



for 
for 



Wsf < T, 



(44) 



which is independent of ^ for ujsi <, T (under-doped re- 
gion). This means that p oc T is expected for wider 
range of filling n, which is consistent with experiments. 
In §V1.B, we calculate p accurately according to the Kubo 
formula. 

At the hot spots, however, A(fc||) deviates from T- 
linear behavior. In fact, A(A:||) at the hot spot is given 
by using as 



A hot — 



3C^ 
47r 



FS 



4wq(a;q + ■kT/2) 



11 



2\v\ 
As a result, 



Ahot oc 



-XQ 



Tr'- (l-(l + ^T/2c^sO~'/'). (45) 



cx VT 

cx VT 



for tJst > T, 
for Wsf ;S T- 



(46) 



Thus, we find that Ahot VT for wider range of filHng 
n, which seems to be realized in Fig. |3| in all the cases. 
This relation does not contradict with the T-linear resis- 
tivity because p is determined mainly by the cold spot 
properties. 

Finally, we show the temperature dependence of the 
anisotropy of A(fc||), r = Acoid/Ahot in Fig. In 
all cases, r becomes smaller in the under-doped region, 
which is consistent with recent ARPES experiments, fs^ j 
However, we see that r depends on the shape of the FS's 
sensitively. The relation r cx VT, which is expected ac- 
cording to eqs. ( ^ ) and (^5|), is satisfied clearly only in 
YBCO. In conclusion, the relation r cx VT is less uni- 
versal than the Curie- Weiss behavior of XQ in HTSC, 
which is reproduced by the FLEX approximation for all 
compounds, (see Fig. |l2|). 



To obtain a^^ and Uxy/H, we solve the following equa- 
tions self-consistently: 



^..^e^E/v (-^) ( |Gk(e)rt^k.(e)Jk.(6) 

-Re{Gl{e)-vl,{e)}), (47) 



x|Gk(e)r^(k,e) 
^^(k, e) = Uka;(e) 

+ {x ^ y), 



dky 



cth th — 

2T 2T 

2 



(48) 



(49) 



xIniVk-q(e - Lu + iS) ■ \G^ie)\' ■ Jq^(e), (50) 
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FIG. 15. The temperature dependence of Acoid/Ahot- 



where /(e) = (exp((e - ^)/r) + l)-i, and G(tJ + i(5) and 
E(cij-|-i(5) are derived from G(w„) and S(ct;„) through the 
numerical analytic continuation. |]53|| The e-integration 
in the above equations are not difficult because its lead- 
ing contribution comes only from |e| <, T. Note that 
|Gk(c^)P = 7rpk(c^)/Ak(c^) and |Gk(c^)|2|ImGk(c^)| « 
^Pk(w)/2A2(a;). 

As for the resistivity, the existence of the second term 
of eq. (p7[), whose derivation will be publised elsewhere, 
has been overlooked in literatures so far. It gives qualita- 
tively an important contrivution in the case of ZkAk ~ 
in where the quasiparticles are not well-defined. 



B. Resistivity and Hall coefficient 

In §V, we find that the vertex correction from the MT 
term gives singular behaviors. In this subsection, we 
obtain the self-consistent solution for Uxx and Uxy/H, 
by using the self-consistent Green function obtained by 
the FLEX approximation. We solve the BS equation for 
Jk^(u;) explicitly by including all the MT terms. Here we 
do not use eqs. (^, ( p^ ) and ( ^ ) because the energy- 
integrations in deriving them have been done under the 
assumption that ZfcAk ^ T. In reality, ZfcAk ^ T is 
realized in the present case. Below, we perform the 
energy-integration seriously by taking account of the 
energy dependence of v\^^ji,{lu) = w^^^ + dKeT,\s_{uo) / dk^j_, 
Ak(a;) = -IniEk(t^ + i*^), Vfe(cj), and Jk^(cj). 



Because Sk(w) satisfies the self-consistency condition, 
v\if_i{uj) includes all the vertex corrections (a)-(c) in Fig. 
^ automatically. Although J\f,{Lo) contains only the (a) 
process of 7kk' (e, e') in the present calculation, the oth- 
ers give only tiny corrections as shown in §IV. In this 
sense, our theory satisfies the condition of the conserv- 
ing approximation well numerically. We did not find any 
difficulty in solving the BS equation ( |50| ) for Jk^(w) nu- 
merically, since the self-consistency condition for Gk('£) 
is satisfied in the FLEX approximation. |2^,^ Figure [iq 
shows the obtained J(fc||) for YBCO on the FS along the 
path shown in Fig.^ Its feature is close to the schematic 
one in Fig ^. We see that Jy is negative around the hot 
spots in this figure. Apparently, such a region is enlarged 
in the case of NCCO. 
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FIG. 16. The obtained J(fc||) together with v{k\\] on the FS. 



Below, we examine the obtained numerical results for 
YBCO, NCCO, and LSCO. The calculated temperature 
dependence of p = ^joxx and i?H — {cTxy/H) ■ (? are 
shown in Fig. |l^ and Fig.nS], respectively. In these fig- 
ures, we also plot = l/cr^ and i?^ = {a^y/H) ■ 

is given by replacing J-k^{uj) with 

„„A ^0 



where cr°^ and a^y 

Vkfj.iuj) in eqs. @ and (|4^). Both a'^^ and cr^^ - cr„„ 
are the result of the conventional Boltzmann approxima- 
tion, where the conservation laws is violated. 
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Resistivity : 

At first, we discuss the T-dependences of the resistivity 
shown in Fig. |l^, where we put e^/h = 1. T — 0.1 
corresponds to ~ 500K if we assume |to| ~ 0.5eV. In ev- 
ery case, both and p show an approximate T-linear 
behavior, reflecting the temperature dependence of Ak 
at the cold spots as shown in Fig. |l3|. |2^,Q They are 
consistent with experiments. In all the cases the rela- 
tion /5 > p° is realized, as is expected from the analysis 
in §V. In LSCO and YBCO, the extrapolated value of 
p'^ at r = from the higher temperature region is zero, 
while that of p seems to take a finite value even in a pure 
system. This behavior of p can be explained by looking 
at eq. because decreases as T increases, reflect- 
ing the decrease of the backward scattering processes at 
higher temperatures. 



The doping dependence of p in YBCO and LSCO is 
very small by using the present set of parameters. Ac- 
cording to Ref. m, however, {p{300K)-p{100K))xlO-'^ 
[flcm] in La2-5Sr5Cu04 is about 4.0, 2.5 and 1.8 for 6— 
0.11, 0.18 and 0.28, respectively. This discrepancy with 
experiments is an important future problem, which will 
be improved by choosing a more appropriate set of pa- 
rameters. Besides, we can not reproduce the doping de- 
pendence of the residual resistivity observed experimen- 
tally because we neglect the impurity effect. We note 
that i?H is more insensitive to parameters than p is be- 
cause i?H is essentially independent of the life-time of the 
quasiparticles. 
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17. Temperature dependence of p and p~ 
Note that p = 1. 
*f2cm in single layer compounds. 'VC 



in all the cases. Note that p = 1.0 in this figure cor- 



FIG 
p > p' 

responds to ~ 4 x 10 
means the vertex corrections for the current. We stress that 
dp/dT increases below T ^ 0.08 in YBCO and LSCO, while 
dp^ /dT does not. This phenomena is caused by the VC, not 
by the psuedo-gap formation in the DOS. 

Hall Coefficient : 



Next, we discuss i?H shown in Fig. At higher tem- 
peratures (T ~ 0.2), where ^ < 1 is satisfied, we see 
that i?H ~ Ryi for all compounds. As the temperature 
decreases, iZn increases for YBCO and LSCO, and de- 
creases for NCCO, and |i?H — -RhI follows the Curie- Weiss 
like behavior in all cases. Moreover, its coefficient in- 



13 



creases rapidly as the filling approaches to n = 1, which is 
consistent with the experimental relation i?H oc |1 — 
These behaviors are consistent with the analysis in §V. 
The obtained filling dependence of i?H is much stronger 
than that of p (or dp/dT). 



Moreover, the sign of i?H in NCCO changes to negative 
below T < 0.08 ~ 400K, which is consistent with exper- 
iments. The Boltzmann approximation can not explain 
this behavior because the shape of the FS is hole-like ev- 
erywhere. In the case of LSCO, the FS around (±7r, 0) or 
(0, ±7r) are convex in the over-doped region. In this case, 
the Hall coefficient can be negative within the Boltzmann 
approximation. Experimentally, i?H in Lai-xSr2:Cu02 
becomes negative and almost temperature independent 
for X >, 0.32, where no superconducting transition occurs 
and the AF fluctuations are much weaker. In the present 
calculation for LSCO at n — 0.65, we find that i?H is 
nearly zero for T = 0.2 ^ 0.02, and Rn « i?^ is real- 
ized. As a result, important features of i?H in LSCO are 
reproduced in the present study for |1 — n| > 0.1. 



Here, we consider the T-dependence of the Hall co- 
efficient given by neg lecting the vertex corrections, 
According to Fig. |ia, decreases moderately in LSCO 
and NCCO, while of YBCO shghtly increases at T < 
0.05. First, we discuss the effect of the T-dependence 
of r = Acoid/^hot on the Hall coefficient. Figure [l5| 
shows that r decreases as the temperature decreases, and 
rocVTis reahzed well for YBCO at n = 0.90. This effect 
makes larger at low temperatures because the effec- 
tive number of carriers, which contributes to the trans- 
port phenomena, decreases around the hot spots. This 
mechanism has been pointed out by several authors to 
explain the enhancement of R^ in YBCO. [HHH] Sec- 
ondly, we discuss the effect of the T-dependence of the 
FS. As shown in Fig. |l^, the curvature of the FS around 
the cold spots decreases as T decreases, which should 
make ji?!^! smaller. In conclusion, through the cancella- 
tion of these two effects, the Hall coefficient given by the 
Boltzmann transport approximation is nearly 1/ne and 
is not enhanced significantly, if we take the temperature 
dependence of the FS into account correctly. 



We stress that i?H in our calculation follows the Curie- 
Weiss law, even if the T-dependence of the FS is taken 
into account. Undoubtedly, this behavior of i?H comes 
from the vertex corrections for the current. In fact, in 
the present calculation, the relation i?H oc XQ given by 
eq. (^ij) seems to be satisfied qualitatively for all the 
sets of parameters. In summary, the vertex correction 
for the current is essential for the Curie- Weiss behavior 
of i?H in high-Tc cuprates observed experimentally. This 
universal behavior of i?H is quite robust in the present 
calculation. 
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FIG. 18. Temperature dependence of Rh and J?h. We see 
that i?H (more precisely Rh — Rh) follows the Curie- Weiss 
type law in all the cases. This universal behavior is ascribed 
to the T-dependence of Here we put e = 1. Note that 
l/\ne\ ~ 1.5 X lO-^'cmVC in HTSC's. 



Hall Angle : 

We also discuss the temperature dependence of the Hall 
angle Oh, which is defined by cot(0H) = (^xx/ {o'xy/H) = 
p/Ru- Figure [ig| shows that cot(0H) is approximately 
proportional to 7^ for 0.02 < T < 0.08. This relation has 
been observed experimentally in various kinds of HTSC's 
for T = 100 ~ 300K. Our theory can explain this 

relation without assuming the non-Fermi liquid ground 



14 



state which possesses two kinds of relaxation rates. |5^] 
Figure |l^ means that i?H follows the Curie- Weiss behav- 
ior, because p is proportional to T. We stress that the 
relation cot(0H) oc is also observed experimentally in 
K-BEDT-TTF compounds psl or in V2O3, (7| both of 



them are also nearly AF- Fermi liquids. 
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FIG. 19. T-dependence of the Hall angle as cot(OH) vs 
in YBCO and LSCO, for T = 0.02 ~ 0.08. (i.e., 
T = 120 ~ 480K.) The thin line represents the cot(eJ![) ob- 
tained within the Boltzmann approximation. 

Here, we discuss the momentum dependence of the fol- 
lowing functions: 



dfcipk(0)yl.(k,0) 



{Ak(0)}^ 



dk 



{A(fc||)} 



(51) 



where fc|Lis the momentum along the FS. As(k, e) is given 
by eq. (|49|), and fcx is the momentum perpendicular to 
the FS. It is clear that <7xy/H oc ^^^dk\\Sxy{k\;^). We 
also define S%y{k\\) by replacing Jk^ with V]^^ in eq. (49), 



which means that a'^y/H oc Jpg dk\\S^y{k\^). 

Figure pO| shows the momentum dependence of Sxy{k\\) 
and 5'^y(fc|t) along the FS shown in Fig. |lj. In both cases 
of YBCO and NCCO, 5*^,^(^11) is positive everywhere. 
Whereas, Sxy{k\\) is positive inside the MBZ and nega- 
tive outside of it, which is consistent with the analysis in 

§v. 

In the case of YBCO, Sxy{kii) takes a maxi- 
mum value on the XY-axis i.e., on the cold spot. 
It takes an enhanced value because the relation 
Sxy{kco\d)/S°y{kcoid) = 1/(1 - a^(fccoid)) oc ^2 is ex- 
pected according to eq. (^. As a result, i?H oc is 
realized. We have also calculated Sxy{k\\) for LSCO, and 
found that its behavior is similar to that for YBCO in 
spite of the incommensurability of x('J'jO). On the other 
hand, in the case of NCCO, Sxy{k\\) takes a maximum 
value on the BZ-boundary which is a cold spot of NCCO. 
It is also enhanced because Sxy{kcoid) / Sxyikcoid) oc 



/3||(fccoid) oc is expected by eq.(^). As a result, 
i?H oc — is realized and i?H becomes negative at low 
temperatures. 
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FIG. 20. The momentum dependence of 5'a;j;(k) and ^^^(k) 
on the FS. 
(see Fig. 



VII. SUMMARY AND DISCUSSIONS 



First, we outline the main results of this paper: We 
have calculated the conductivity axx a-nd the Hall con- 
ductivity Gxy/H in the single-band Hubbard model based 
on the Fermi liquid theory. We have calculated the total 
current Jk including the vertex corrections in the conser- 
vation approximation. In nearly AF Fermi liquids, the 
Bethe-Salpeter equation (|l|) for Jk can be simplified to 
eq. ( ^5|) approximately. The obtained Jk shows non- 
trivial critical behaviors as seen in Fig. ^, which is the 
natural consequence of the strong backward scatterings. 
In conclusion, Ra oc is realized in HTSC's through the 
anomaly of Jk. This mechanism has not been pointed out 
previously. 

We also have done the numerical calculations by using 
the FLEX approximation. We can reproduce character- 
istic features of the spin fluctuations for YBCO, NCCO 
and LSCO, by using the appropriate set of parameters. 
In each cases, dx2_y2 superconductivity is realized at 
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Tc = 50 ~ lOOK. Next, we have determined Jk by solv- 
ing eq.(|l8|) numerically, and calculated both p and i?H 
for various filling numbers. As shown in Fig. the 
overall features of i?H in each compounds are reproduced 
quite well. Especially, both the relations i?H oc 1/T and 
p oc T are obtained at the same time. We have found 
that i?H < is realized in NCCO because the cold spots 
in NCCO locate around the BZ-boundaries, which may 
be verified by ARPES experiments. 

The vertex corrections mentioned above are not in- 
cluded in the Boltzmann approximation. We have con- 
firmed that the Hall coefficient given by the Boltzmann 
approximation remains of order 0{l/ne) if we take 
the T-dependence of the FS into account correctly, (see 
Fig.[l^) . Moreover, remains positive because the FS 
is hole-like everywhere. In conclusion, the anomalous be- 
haviors of i?H in HTSC is reproducible only if the vertex 
corrections for the current are taken into account. 

Here, we discuss the validity of the relation Rn oc 
In a conserving approximation (including the FLEX ap- 
proximation) the interaction Vk{uj) which gives ImEk(w) 
also determine the MT-type vertex corrections for Jk, as 
shown by eqs. (|3ll) and (p3|). This condition leads to 
ttk ~ 1 in eq. (pSj), which strongly suggests that the 
relation i?H cx; f^vill be valid beyond the FLEX ap- 
proximation. Now we assume that it is valid near the 
half-filling case (n w 1). Then, the experimental rela- 
tion max{|_RH|} oc 1/|1 — n| can be understood because 
max{^^} cx 1/|1 — ri| is observed experimentally near the 
half-filling. Next, we consider the Hall coefficient be- 
low the pseudo spin-gap temperature, T*, where i?H in 
YBCO decreases as T decreases experimentally. It is also 
understandable because ^ slightly decreases below T* ex- 
perimentally. 

We also find that the resistivity is slightly enhanced 
by the vertex corrections for the current. Moreover, we 
obtain the expression of the incoherent conductivity (Tine, 
which is given by the second term of eq.(^). In the case 
of ZkAk ~ T, which is satisfied in HTSC's, Oinc (> 0) can 



be the same order of the first term of eq.(47) at higher 
temperatures. 

Unfortunately, the FLEX approximation becomes in- 
sufficient near the Mott-insulating state. By this reason, 
we did not apply the present method for 0.9 < n < 1.1. 
Experimentally, both |_Rh| and dp/dT for 0.9 <n< 1.1 
increases rapidly as n ^ 1. The FLEX approximation is 
also inappropriate for the study of electronic states be- 
low T* , which is one of the important future problems on 
HTSC. Recently, Ref. @ calculated some vertex parts 
of the self-energy. As an alternative possibility, the pre- 
formed pairs may be formed for T < T* . This scenario 
has been intensively studied recently. p9| , |60[ 

In Ref. [0 , it is shown that a similar numerical study 
based on the AF spin-fluctuation model also leads to the 
enhancement of i?H, by using the set of spin fluctuation 
parameters consistent with experiments. Although the 
conserving laws are not satisfied exactly in that study. 



it confirms the importance of the vertex corrections for 
the current. It also indicates that the numerical results 
presented in this paper should not be taken as an artifact 
specific to the FLEX approximation. 

Now we would like to discuss i?H in heavy Fermion 
(HE) compounds. In various paramagnetic compounds, 
i?H shows a drastic temperature dependence and takes 
an enhanced value. ||6l|-|63[| At low temperatures, the re- 
lation i?H = c • is observed in many compounds, and c 
is always positive. It is explained in terms of the anoma- 
lous Hall effect (AHE), which originates from the local- 
ized /-orbital angular momentum, and its enhancement 
factor is given by xo = Xq=o(0). p^ , |65[ On the other 
hand, in some HF compounds with AF ground state, 
the relation i?H — c ■ p^ (c > 0) is not satisfled and 
the sign of dRu/dT at T > Tn changes depending on 
compounds. ||6l|-|63[| Similar behavior is also observed in 
non- Fermi liquid HF compounds, which is near the AF 
quantum phase boundary, e.g., Ce(Nii_2;Pda;)2Ge2. |6^] 
We stress that the normal Hall effect can exceed the AHE 
and i?H oc XQ will be realized in these nearly antiferro- 
magnetic HF compounds, where xq S> Xo is realized. 
We note that XQ oc (T — Tn)^'^/^ in three dimension is 
obtained by the SCR theory. 

We also comment on the k-BEDT-TTF organic super- 
conductors. The recent studies by the FLEX approxima- 
tion reveal that the large AF fluctuations are the origin 
of the d-wave superconductivity. Recent experi- 

ments show that the _Rh of this system increases as T 
decreases |]67| , |68| , which can be reproduced well by us- 
ing the analysis of this paper. ||6^ Also, the relation 
i?H oc 1/T is observed in the superconducting ladder 
compound Sri4_2:Caa;Cu2404i. |70|, whose characteris- 
tic electronic properties are well explained by the FLEX 
approximation. pSj] Moreover, i?H in V2O3 shows the sin- 



gular T-dependence near the AF phase boundary. [171 72 
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APPENDIX A: ANALYTIC CONTINUATION 

FOR r(eiv,ejv') 

In this appendix, we derive the irreducible vertex 
T]ik'{e,e.') which are the kernel of the BS equation, eq 
(p^). For this purpose, we perform the analytic contin- 
uation for r(e„, e„'; cj;), where ui — 2'kTI (I > 0) is the 
external frequency. The irreducible vertices consistent 
with the FLEX approximation are given by (p5|)-(p7|). 
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According to eq. (12) in ref. ||27[|, 



+cth^(r™(e,6')-rLX-(6,6')) 

-tii^(rLV(e,0-rS:'(e,e')). (Ai) 

Here, r"(e,e'), r™(e,e') and riV(e,e') are given by the 
analytic continuations of r(e„, e„/; wj) for regions II, 111 
and IV in the complex (e, e')-plane shown in Fig. pl| , 
respectively. 
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FIG. 21. r(e, e'; ui) (uj > 0) is an analytic function inside of 
each (II,III,IV)-region, and has the cuts on each lines. 

Next, we take the limit icui +0. For eq. (p5|), we get 



r^^^}\e,e')^V^ik^k',e'-e), (A2) 



l/^(k-k',e'-6). 



(A3) 



where A{R) represents the advanced (retarded) function. 
Taking account of the relation Im{e' — e} > in the II- 
region and Im{e' — e} < in the llLIV-region, we get for 
eq. (p6|) as 



th'^lmG^+Ju; + e)GS,+Jcj + e') 



- th 

-C. 



2T 



(A4) 



th^^:^ImG^+q(c. + e)Gi^,+q(c. + e') 

+ th^^G^+q(c^ + e)ImG^,+q(c. + e') 
+G, (A5) 
where G is a real function, and Wq{uj) is given by 



+ -l/^|C/X^(a;)-l| -U\ 



(A6) 



In the similar way, we can obtain the expressions of 
r|;^"(e,e') = Tl^l}'\e,e') and tI^I]"" {e,e') through the 



analytic continuation of eq. (|2j), by taking account of 
the relations lm{e' + e + w} > in the 11, Ill-region and 
Im{e' -I- e -|- w} < in the IV-region. 

By inserting the above equations into eq. (|Al|), 
Tk,k' (£„,£„' ;w/) is given by T(") + T(^) + T(^\ They 
are derived as 



cth- 



2T 



thil 

2T 



x2i ImVk'-k(e' - e + i(5), 



(A7) 



cth- 



2T 



X —ITT th th 

' \ 2T 2T 

xpk+q(w + e)/3k'+q(t^ + e'): (A8) 
e' + e 



cth- 



2r 



x(-i7r) th-::^ +th 



2T 2T 
X/Ok+q(w + e)/9k'-q(-w + e')> (A9) 



Note that T\^^''^\^i^') is pure imaginary. 

APPENDIX B: CORRECTIONS FROM THE 
AL-TERMS FOR Jk^ 

In this appendix, we study the contributions from the 
AL-terms to the total current Jk- For this purpose, we 
solve the BS equation for Jk(w) including both the MT 
process and AL precesses, and compare the obtained re- 
sults with those in §VI B. The exact Bethe-Salpeter equa- 
tion is, 

a,6,c 

JM=v^,,{uj) + Y,^JL^^). (Bi) 

r 

/OO 7 
^^T'it^^^^)-\Gq{e)\^-Jq,{e), (B2) 



where r = a,b,c and ^^q(w,e) are given by eqs. (A7)- 

(|A9|). Note that ^'q(w,e) are pure imaginary. 

For simplicity of the numerical calculation, we put all 
the energy variables in pit,'{io), Ai^(uj) and Jit,fj,{uj) as zero 
for AJ^^{uj) and AJ^^{cj). Strict justification of this sim- 
plification is difficult, although it may be sufficient for a 
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rough estimation of the magnitude of the AL terms. We 
represent the sohition of (Bl) as J^/J^^^i^)- 

Figure shows the calculated resistivity and the 
Hall coefficient for YBCO (n = 0.90) derived from 
J^'^"'"^^(u;), together with those given in §VI B. We see 
that the AL terms give only small corrections to p and 
i?H. 
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FIG. 22. 

The obtained T-dependence of the resistivity and the 
Hall coefficient. MT and MT+AL are given by Jk;^(w) 
derived from eq. (|50|) and (Bl), respectively. 
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